Scattering of p/x muonic atoms in solid hydrogen 
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We present the results of experimental and theoretical study of the scattering of low energy pfi 
atoms in solid hydrogen cooled to 3 K. Strong effects resulting from the solid state interactions have 
been observed in the TRIUMF experiment E742 where muons were stopped in thin frozen layers of 
hydrogen. The resulting emission of low energy p^t atoms from the hydrogen layer into the adjacent 
vacuum was much higher than that predicted by calculations which ignored the solid nature of the 
hydrogen. New differential scattering cross sections have been calculated for the collisions of pfj, 
atoms on solid hydrogen to account for its quantum crystalline nature. Analysis of the experimental 
data performed using such cross sections shows the important role of the coherent scattering in 
Pfj, atom diffusion. For p/i energies lower than the Bragg cutoff limit (~ 2 meV) the elastic Bragg 
scattering vanishes which makes the total scattering cross section fall by several orders of magnitude, 
and thus the hydrogen target becomes transparent allowing the emission of cold p/i atoms to occur. 

PACS numbers: 36.10.Dr, 39.10.+J, 61.18.Bn, 34.50.-s 
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I. INTRODUCTION 

Negative muons stopping in hydrogen can form muonic 
hydrogen (p/i) atoms. Although created in excited 
states, such atoms cascade to the ground state quickly 
(10~ 12 s) where their kinetic energy is of the order of 
several eV, much higher than thermal equilibrium en- 
ergies. The muonic hydrogen atom is about 200 times 
smaller (m^/im, scaling) than the size of ordinary elec- 
tronic hydrogen. The small neutral atom can easily dif- 
fuse through the surrounding medium undergoing differ- 
ent types of interaction including elastic and inelastic 
scattering. Scattering of fast p/x atoms in hydrogen is 
governed by a large cross section (a s > 10" 19 cm 2 ) which 
is quite effective at slowing them down. 

Only a few experiments have examined the scatter- 
ing of muonic atoms on nuclei and molecules directly, 
although it is an important process in most muon 
physics phenomena such as muon catalyzed nuclear fu- 
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sion (/iCF) or muon nuclear capture by protons (see re- 
views P, 0, IH 0) The first experiments of p/i + H 2 
scattering were performed in gaseous hydrogen and used 
the traditional diffusion method 0, IE 0, 13 • 

Much more has been done to study the scattering the- 
oretically. Many calculations of cross sections for scat- 
tering on bare nuclei, atoms, and hydrogen molecules 
have been made, however, solid state effects were not 
considered. The high accuracy calculations of the total 
cross sections for low energy scattering (e co ii < 50 eV) 
for p/i and other muonic atoms on bare hydrogen nuclei 
(called the "nuclear" cross sections) were done in Ref. [jj 
by solving the Coulomb three-body scattering problem 
using the adiabatic multichannel approach. Differential 
cross sections for that case were calculated in Ref. [lfj 
using phase shift values from Ref. 0. For collision ener- 
gies lower than about 0.1 — 1 eV it is necessary to account 
for both electron screening and the molecular structure 
of the target. Total and differential cross sections for this 
case (called the "gas" cross sections) are given in Ref. ^l| 
and in Ref. 0] , respectively. Another possible approach 
to including the molecular effects for epithermal energies 
uses the Sachs-Teller tensor-of-mass model and can be 
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found in Refs. dill. 

The scattering experiment results given in Refs. 0, ||| 
0, 0] were sometimes inconsistent not only among them- 
selves but with theory as well (see ^3 f° r a review) . The 
latest and most advanced measurements in gaseous hy- 
drogen were performed at PSI where the cross sec- 
tions for p/i scattering on H2 molecules [TT| were used 
for the analysis of the experimental data. Those mea- 
surements were not in full agreement with the theory. 
On the other hand, the d/x + D2 scatte ring measurements 
performed by the same collaboration [16j were in agree- 
ment. 

Until now, no experimental studies concerning //-atom 
scattering in solid hydrogen have been performed. Such 
experiments are complicated to analyze because the re- 
sults for the cross section are not directly obtained but 
are only deduced by their effect on other results, such as 
time distributions or yield intensities, which themselves 
are often obscured by other background processes. 

The development at TRIUMF of the multila yer thin 
frozen hydrogen film targets [H H |U M, IH H3, 
which produce muonic atom beams emitted into vacuum, 
permitted the cross sections to be probed in another way. 
We have studied several isolated muon induced processes 
using a time-of-flight (TOF) method permitted by the 
frozen target geometry pflj . In particular it was used in 
TRIUMF experiment E742 for the cross section study of 
d/i + H 2 , t/i + H 2 scattering and the Ramsauer-Townsend 
effect (RT) which is seen for these systems at collision en- 
ergies between 2 — 10 eV p3.l25ll26| . During those mea- 
surements a strong emission of low energy p/i atoms from 
the hydrogen layers into adjacent vacuum was observed. 
The yield was much higher than expected based on cal- 
culations which ignored the solid nature of the hydrogen 
target. Additional experimental studies H3,|2^| and new 
theoretical calculations of "solid" cross sections [2^ have 
been performed in order to clarify and explain the p/i 
emission behavior. 

This paper summarizes our findings. In section II the 
theoretical background and new calculations of scattering 
cross sections in solid hydrogen are described. The ex- 
perimental apparatus and the measurement method are 
given in section III. Section IV presents the measurement 
results and their analysis, whereas section V contains the 
discussion of the results and some concluding remarks. 
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FIG. 1: Energy levels for the p/i + p system, where £1 and 
£2 represent the collision energies for the singlet (F = 0) and 
triplet (F = 1) states, respectively. The dashed line repre- 
sents the total energy of the system. Four transitions with 
cross sections i,k — 1,2 are possible due to the hyperfine 
splitting of the energy levels. For the scattering of p/i in the 
singlet state only the elastic scattering an is possible if the 
collision energy £1 is below Ehfs = 0.182 eV. 



as a quantum mechanical Coulomb three-body problem. 
Figure ^ shows the diagram of the two isolated states of 
the p/i + p system with the possible transitions. 

Most of the following calculations were performed in 
the adiabatic representation which results from expand- 
ing the wave function of a three-body system over a com- 
plete set of solutions of the two-center Coulomb prob- 
lem |33, especially in the two-level appr oximation |32] | 
with further modifications 02 HE HH • The progress 
in the effective potential calculation for the two-center 
problem permits the multichannel scattering equations 
to be solved even when there are a large number of closed 
channels [3^ . l38| . Accurate calculations of the total cross 
sections performed in this multichannel approach includ- 
ing p/i + p scattering are presented in Ref. 0. Reac- 
tance T-matrices and phase shifts, also given in Ref. [jj, 
for different values of the total orbital angular momen- 
tum and the total spin of the three-particle system have 
been used to calculate the differential cross sections [lOj . 
Cross sections (total and differential) for collisions with 
energies less than w 1 eV, where both electron screening 
and molecular binding are important ("gas" cross sec- 
tions), are given in Refs. [HEl. The screening effect is 
described there in terms of the effective screening poten- 
tial, and the Fermi pseudopotential method was applied 
to model the chemical binding. 



II. SCATTERING CROSS SECTIONS 
A. Scattering on bare nuclei, atoms and molecules 

Cross section calculations for the IS state p/i atom 
scattering on bare protons 

PM(F)+ P ^pm(F')+P, (1) 

(where F, F' are the initial and final muonic atom spins) 
were begun by Gershtein [3(j. He treated the process 



B. Cross sections for solid hydrogen 

The TRIUMF experiments E El IH HI |H H| 

have stimulated theoretical studies of /x-atom scatter- 
ing in solids. Solid hydrogen at zero pressure is a quan- 
tum molecular crystal, which is characterized by a large 
amplitude of zero-point vibrations of the molecules. At 
3 K, the vibration amplitude is approximately 18% of 
the nearest neighbor distance for the H2 molecule, and 
15% for the D 2 molecule [3^. Nevertheless, experiments 
show that quantum solids display typical crystal struc- 
tures and that common crystal characteristics, such as 
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the density of the vibrational states, are well defined. 
This proves that the molecular motion is correlated in 
such a manner that the crystalline structure is not de- 
stroyed. However, theoretical methods developed for a 
classical-crystal description, encounter certain problems 
when applied in the case of quantum crystals. Namely, 
the interaction potential between the hydrogen molecules 
has a highly repulsive anharmonic core and thus the stan- 
dard lattice dynamics leads to imaginary vibration fre- 
quencies. Nevertheless, the standard dynamics can be 
used, after a renormalization of the interaction potential 
by accounting for the short-range correlations between 
neighboring-molecule movement. 

Solid hydrogen can exist in both efficient-packing 
structures: face-centered cubic (fee) and/or hexagonal 
close packed (hep). Since the TRIUMF targets were 
formed by rapidly freezing hydrogen gas on a gold foil 
at 3 K and zero pressure, the solid layer has a polycrys- 
talline fee structure |39j | . There is also experimental evi- 
dence that thin hydrogen films formed on fee metal (e.g., 
gold or silver) will retain that same structure ^3] ■ Since 
the fee and hep crystals are very similar (e.g., the molar 
volumes are almost the same and the first three shells of 
neighbors of any fixed molecule are identical in both the 
structures) , the fee cross sections are also a good approx- 
imation to the hep case. The hydrogen was purified by a 
palladium filter at 600 K immediately prior to freezing, 
so the resulting solid target had a statistical distribution 
(1:3) of molecular rotational states K = and K = 1. 
Such a mixture of rotational states is often called "nor- 
mal" hydrogen, nH2. The relevant lattice constant for 
the fee structure at zero pressure is 0.5338 nm |3Sl | - 

A method to calculate the scattering cross sections of 
muonic hydrogen atoms in solid hydrogen ("solid" cross 
sections) based on Van Hove's approach and using phase 
shifts for muonic atom scattering on bare nuclei 0, 
has been proposed by Adamczak (see [2|j for details). 
The calculated differential cross sections include incoher- 
ent and coherent effects. The impinging muonic atom 
can induce inelastic reactions, both in a single molecule 
(rovibrational transitions, or spin-flip) and in the whole 
target (excitations or deexcitations of the lattice vibra- 
tional states). The latter are usually interpreted as cre- 
ation or annihilation of phonons. It is possible to create 
or annihilate one or more phonons in coherent or incoher- 
ent processes, but in practice, annihilation processes are 
strongly suppressed in a 3 K target because few phonons 
exist at low temperatures. 

Figure [3 presents the calculated total cross sections 
for p/n scattering on 3-K solid fee hydrogen for differ- 
ent initial and final spin states F of the p/i atom. Also 
are shown some details of the total cross section for p/i 
atom scattering from the ground spin state F = 0. For 
the sake of comparison, the doubled cross section cr™ c of 
p^(F = 0) + p nuclear scattering is plotted. The cross 
sections are given for a single bound molecule. At ener- 
gies greater than roughly 1 eV both the solid state and 
molecular binding effects are very small and therefore the 



o 



S3 
O 

- rH 

o 

CD 

t/5 

VI 
t/5 
O 

& 10 



10 
10 
10 
10 

1 



— 1 — 1 — 1 — 1 — 

r^-^a 21 p|j.(F)+nH 2 . T=3 K 


: a 22 
\ a n 






a i2 

■ ■■■i 


: 2xa™ c 


)N loss 

; o„ 

sain 

1 ' ■ • ■■ •. a 
•„■:.-.. ■: „ 1 1 

i . V r , . 



10 

10 4 10 3 10 2 10 1 1 10 

p(0, energy (eV) 

FIG. 2: Total cross sections for p/i(F) scattering in 3-K poly- 
crystalline nH2 with the fee structure, for different values of 
the initial and final muonic atom spin, F. The dotted line rep- 
resents the phonon-annihilation fraction of an that results 
in p/i energy gain; the sum of contributions from phonon cre- 
ation and rovibrational excitations to an, which lead to p/i 
energy loss, is denoted by dash-dotted line. The doubled nu- 
clear scattering cross section a™ for p/i(F = 0) + p is shown 
for comparison (dashed line). Note the Bragg cutoff energy, 
E_b, at e ps 2 meV for an. 



cross section for a real hydrogen target is well described 
by the nuclear cross section. 

In Fig. |2J there is an important difference between the 
singlet (iTii, F = 0) and triplet (022, F = 1) state scat- 
tering. For the singlet, only the state J = | of the total 
spin of the p/i + p system is possible. As a result, the 
scattering in nH2 is almost fully coherent and thus inter- 
ference effects determine the behavior of the singlet cross 
section at the lowest energies. Below the Bragg cutoff 
energy, E# ~ 2 meV, elastic and phonon-creation coher- 
ent scattering is impossible and the total cross section 
is determined by the weak incoherent processes, which 
gives rapid falloff of an. Coherent phonon annihilation 
is present below Eg, but its magnitude is very small at 
3 K. The rotational deexcitation K = 1 — ► K' = of 
an H2 molecule gives no contribution to the cross section 
since this transition is strictly forbidden for F = 0. 

Scattering of the p//(F = 1) atom is possible in the two 
total-spin states: J = h and J = |. The nuclear am- 
plitudes for these two processes are very different and 
therefore averaging the molecular scattering amplitude 
over spins leads to a strong incoherent component. As a 
result, the cross section 022 m solid hydrogen has a large 
magnitude at e — > 0, though a small falloff of its value 
at Eb is still present. Significant contribution to ayi at 
lowest energies comes also from the rotational deexcita- 
tion K = 1 — > K' = 0, which is possible (for F = 1) due 
to the exchange of the muon between the protons during 
the collision process. 

In Fig. [21 are shown contributions to <7n from different 
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processes. The energy region 2 — 10 meV is dominated 
by the strong elastic Bragg scattering. Phonon annihila- 
tion, denoted by the label cxf J 111 , is the only mechanism 
of p/i acceleration. Weak incoherent elastic scattering 
is most important below E#. Slowing of p/i is possi- 
ble through the lattice excitations and then, at sufficient 
incident energy, through subsequent rotational and vi- 
brational excitations. The rovibrational transitions may 
take place with simultaneous one or multiphonon cre- 
ation. The curve which shows the sum of contributions 
from all these processes is labeled by crjf 8 . At e > ujjj, 
the inelastic processes are most important and the cross 
section for the solid (per single molecule) approaches the 
corresponding one for a free H2 molecule. 
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FIG. 4: Slowing down of p/i atoms in a solid hydrogen layer 
of thickness 3.4 mg • cm -2 (1000 Torr-^, as defined in Sec- 
tion [^^XJ. The top plot a) shows the time and energy after 
each scattering event during the whole slowing down process. 
There are 2376 emission events and each p/i atom undergoes 
on the average 40 collisions before emission take place. The 
bottom plot b) shows the time and energy for the p/i atoms 
which have been emitted from the hydrogen layer. The "solid" 
cross sections were used in the simulations. 



FIG. 3: Total cross section for Bragg scattering of p/i(F = 0) 
atoms in 3-K polycrystalline nHb with the fee (solid line) and 
hep (dotted line) structure. 

Figurc[3]illustrates the small differences between Bragg 
scattering of p/i in the fee and hep polycrystalline nH2 . 
The Bragg cutoff energy is slightly (« 0.2 meV) smaller 
in the hep target. Different Bragg peak patterns in the 
total cross sections are distinct only below a few meV. 
The magnitudes of the cross sections are similar in the 
two lattices and, therefore, the theoretical estimation of 
cold-p/i emission, obtained in this work for an fee target, 
is also valid in the hep case. 



C. Slowing down of p/i atoms in solid hydrogen 

The slowing down of p/i in solid hydrogen has been sim- 
ulated by the Monte Carlo method using the new cross 
sections. The more important characteristics are shown 
in the following figures. 

The simulations, performed with the Monte Carlo |42^ . 
represent a real experimental situation where muons were 
stopped and formed p/i atoms in a solid protium target of 
thickness 3.4 mg-cm" 2 (an experiment labeled later as #3 
in Table llV|) . For this presentation only the histories of 



p/i resulting in the upstream emission from the hydrogen 
layer have been chosen. 

Scatterplots on Figs.0]and|3]illustrate the time depen- 
dence of the p/i slowing down (shown are 2376 histories 
ending by the emission from the hydrogen layer). Each 
point represents the p/i energy after a scattering event at 
a given time. Figure^,) shows a sampling of all events 
during slowing down, whereas Fig. 0Jd) shows only the fi- 
nal coordinates when p/i emission has occurred. One sees 
that the slowing down process is very fast and that after 
approximately 10 ns the ~meV energy region is reached. 
Further decelerations of the p/i are then slower since the 
responsible inelastic cross sections become lower. This 
transient region extends to 100 — 200 ns when equilib- 
rium energy is reached and p/i atom diffusion in hydro- 
gen takes place. The equilibrium energy is established 
near the Bragg cutoff limit where both the phonon cre- 
ation and annihilation components of an become equal 
(see Fig. HJ. 

Figure[S]represents contributions to slowing down from 
separate processes for p/i in singlet and triplet states. In 
any experiment both p/i atomic spin states will be ini- 
tially populated, however, the downward spin-flip is so 
fast in the solid target that, after 0.1 ns, practically all 
p/i atoms are in the ground spin state. Therefore, further 
slowing down is governed by the cross section an (Fig. [2J) . 
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FIG. 5: Details of the scatterplot from Fig. namely the 
spin-flip contribution (top) and the non-Bragg scattering pro- 
cesses in the slowing down of p/i atoms (bottom). Note that 
downward spin-flip (F = 1 — > F = 0) increases the energy by 
£hfs = 0.182 eV and the upward spin-flip decreases the en- 
ergy by that amount. Elastic coherent scattering (not shown 
in the figure) is a dominant process for energies near the Bragg 



cutoff energy limit (r 
changed. 



2 • 10 eV) and leaves the energy un- 



Efficient slowing down finishes after about 10 ns and sub- 
sequent p/i diffusion is determined by elastic Bragg scat- 
tering and inelastic phonon scattering. The competition 
between those two processes is shown in Fig. 

The elastic Bragg scattering does not change the p/i ki- 
netic energy in a solid hydrogen (contrary to a gas, where 
the elastic scattering is an effective deceleration process 
because the p/i atom can always transfer a part of its en- 
ergy to a recoiling free H 2 molecule) . Only inelastic scat- 
tering can cause p/i deceleration (or acceleration from the 
phonon annihilation process) but is a weak contribution 
at low energies and at low temperatures. Therefore, p/i 
atoms spend a relatively long time in the diffusion stage 
before reaching the Bragg cutoff energy. In the case of a 
solid H2, rapid falloff of the phonon creation cross section 
at energy Eb makes the p/i thcrmalisation less deep than 
in a gaseous hydrogen. Indeed, the equilibrium energy 
defined by the intersection of ofi™ and a l °f s at ~ 2 meV 
(see Fig. |5J) is still higher than thermal equilibrium in a 
3-K gaseous H 2 . 

The strong increase of the p/i(F = 0) atom mean free 
path due to the sharp decrease of the cross section be- 
low the Bragg cutoff energy is shown in Fig. Such 
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FIG. 6: Competition between Bragg scattering a) and non- 
Bragg scattering b) in p/i collisions in solid hydrogen for low 
energies. Note the logarithmic scales. 
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FIG. 7: Sampled values of the p/i atom mean free path be- 
tween consecutive collisions versus the collision energy. A 
strong increase of the mean free path is seen below the Bragg 
cutoff energy. 



behavior leads to an enhanced emission of cold p/i's from 
the thin solid hydrogen layers. We note that a similar 
phenomenon is used in neutron physics to extract cold 
neutrons from beams produced in nuclear reactors (poly- 
crystalline filters, see e.g., [431]). 



III. DESCRIPTION OF THE EXPERIMENT 

A. The apparatus 

The experiments studying p/i scattering in solid hydro- 
gen were performed at the M20B muon channel at TRI- 
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FIG. 8: The apparatus layout showing the muon entrance 
scintillator (1), the upstream (2) and downstream (3) gold 
foils (themselves inside the cryostat). The surrounding detec- 
tors were the germanium detectors Gl (4) and G2 (5), the 
Nal (6), and the three pairs of electron counters (7). The 
drawing is not strictly to scale. 



UMF. The layout of the apparatus is shown schemati- 
cally in Fig. |H1 Gaseous hydrogen (or neon) was sprayed, 
using a special diffusion system, onto the 51 /im thick 
gold foil, maintained at 3 K, where it froze creating the 
thin solid films which could be maintained in high vac- 
uum. The diffuser was inserted from below and could 
be used to deposit gas on either of the the two gold foils 
separately. The thickness of the film was controlled by 
adjusting the amount of gas injected. Multi-layered tar- 
gets could be made in which the frozen material consisted 
of uniform layers, each made from different hydrogen iso- 
topes or other gases such as neon. A versatile gas han- 
dling system allowed mixtures of different hydrogen iso- 
topes to be prepared with high precision. The frozen film 
deposition uniformity has been measured independently 
via energy loss of alpha particles The amount of 

gas injected into the system was conveniently measured 
in units of Torr • £, where the conversion factor between 
Torr • I and fig ■ cm~ 2 has been determined |44| and is on 
average 3.4 /ig • cm~ 2 per Torr • i for H2. Details of the 
target construction and working procedure are given in 
Refs. ^1,0] • Details of the data acquisition electronics 
can be found in Ref . [2(j ■ 

Incident muons of momentum 26.70 (or 26.25) MeV/c 
(Ap/p = 0.07 FWHM) were detected by a 127/^m scintil- 
lator (1) before traversing a 25 /im stainless steel vacuum 
isolation window. The muons continued to lose energy 



while passing through the cryostat 70 K heat shield to 
eventually stop either in the gold target support foil (2) 
or in the « 400 — 800 /xm thick solid hydrogen target 
where they finally formed muonic atoms. The hydrogen 
target frozen on foil (2), which was placed perpendicu- 
larly to the muon beam axis, was called the upstream 
target (US), and was made of pure protium or of pro- 
tium with a small admixture of deuterium (or tritium) , 
depending on the experiment (see Figs. l9l and fTU|l . 

In the case of a pure protium US target, a thin neon 
layer was additionally deposited on top, as represented in 
Fig. El In other cases, when deuterium or tritium were 
present in the US target, we used an additional down- 
stream target (DS) frozen on a second gold foil placed 
parallel to the first foil but 17.9 mm further along the 
beam axis ((3) in Fig. |HJ. Such an arrangement is pre- 
sented in Fig. where a thin layer of neon is shown 
sandwiched between a layer of pure protium and the 
DS gold foil, and was used for the TOF studies of the 
Ramsauer-Townsend effect. A low muon beam momen- 
tum was chosen to minimize the number of \i stops in the 
DS protium layer. 

Neon was used to detect the scattered muonic atoms 
which left the hydrogen layer and subsequently trans- 
ferred the muon to the neon. The resulting emission of 
207 keV x rays from the 2p-ls /xNe transition was ob- 
served by two ^100 cm 3 germanium crystals (Gl and G2, 
Fig.© with a time resolution of 10- 12 ns (FWHM). The 
Gl detector was used during the whole experiment, dur- 
ing both the deuterium and the tritium measurements. 
However, there were two physically different G2 detec- 
tors, one for each of the deuterium and tritium mea- 
surements. The plastic scintillators ((7) in Fig. were 
located around the target to detect the muon decay elec- 
trons. The Nal detector was used to study /zCF in hy- 
drogen and deuterium mixtures [2^ 0] . 



B. The method 

Following muon capture and p[i formation, the p/x 
atoms slow down and diffuse in the hydrogen layer, with 
some significant fraction of the p/z atoms escaping the 
layer. Analyzing the emission yield and the time distri- 
bution of the escaped p/x's gives information about the 
scattering cross sections. The essential part of the analy- 
sis is the comparison of the experimental yields and time 
distributions with the ones calculated by Monte Carlo. 
The measurements were performed in two different ways; 
either by using the single pure protium target covered 
with a Ne layer (Fig. EJ) or by using the parallel target 
scheme (Fig. ITT))) . 

The first method (Fig.EJ) proved better for studying p/i 
scattering for two reasons: (1) high stopping rate in the 
US target assured high pfx emission statistics, and (2) the 
time spectrum was "clean" in that it did not contain the 
overlapping RT part. Table [TJ summarizes the different 
measurements performed for the study of the RT effect 
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TABLE I: Different measurements performed for the RT and p/i diffusion studies. DEm (Deuterium Emission) stands for a 
layer of 1500Torr -l(H 2 + 0.05%D 2 ) covered with 500Torr-f H 2 . TEm (Tritium Emission) — 2000Torr -£(H 2 + 0.12%T 2 ). sTEm 
(Small Tritium Emission) — 1000 Torr ■ £(H 2 + 0.12%T 2 ). PP(Pure Protium) — 2000 Torr • I H 2 . sPP(Small Pure Protium) - 
1000 Torr ■ I H 2 . GMU — Good Muons: i.e., events when only one muon entered the apparatus (no pileup). Conversion factor 
(for hydrogen): 1 Torr • I corresponds to 3.4 /ig • cm -2 for H 2 . 



Label 


Experimental 


Beam 


US Hydrogen 


US Neon 


DS Protium 


DS Neon 


GMU 




Purpose 


MeV/c 


Torr • £ 


Torr • £ 


Torr • £ 


Torr • t 


xlO 6 


Dl 


RT 


26.70 


DEm 


— 


— 


100 


326.9 


D2 


RT 


26.70 


DEm 


— 


— 


50 


183.3 


D3 


RT, diff 


26.70 


DEm 


— 


300 


50 


521.8 


D4 


RT, diff a 


26.70 


DEm 


— 


600 


50 


433.2 


D5 


diff 


26.70 


DEm 


100 






96.6 


D6 


diff 


26.70 


DEm 


50 






136.9 


D7 


diff 


26.70 


PP 




300 


50 


149.4 


Tl 


RT 


26.25 


TEm 






30 


113.5 


T2 


RT 


26.25 


TEm 






50 


174.2 


T3 


RT, diff" 


26.25 


TEm 




350 


50 


405.3 


T4 


RT, diff 


26.25 


sTEm 




500 


50 


147.1 


T5 


diff 


26.25 


sPP 


10 






199.3 


T6 


diff 


26.25 


sPP 


20 






195.8 



a D4 and T3 are not useful for the pfj, diffusion analysis due to the 
strong overlap between RT and diffusion parts of the time spectra. 



and p/i + p scattering. 

In the second case (Fig. 1100 . the US target was com- 
posed of protium with a small concentration of D2 (or 
T2), which served as a source of energetic d/i (t/i) atoms 
emitted from the layer into the adjacent vacuum with 
a mean energy of 3.5 eV (9 eV in the case of H 2 /T 2 
mixture) as a result of the RT effect. Due to the low 
deuterium (tritium) concentration, p/x's were predomi- 
nantly formed as a result of muon stops in the US layer. 
A small fraction of them survived their evolution in the 
US hydrogen (no muon transfer to the heavier hydro- 
gen isotope, no pp/i formation, or muon decay) and left 
the solid layer after multiple scattering. However, they 
were very slow (e « meV) and could not reach the down- 
stream target before the muon decayed. In most cases 
the muon was transfered from p/i to a deuterium (tri- 




Ne X rays 



tium) atom. At formation, the muonic deuterium and 
tritium atoms had a relatively high kinetic energy, about 
45 eV, which they subsequently lost in elastic collisions, 
mainly with protium, until the energy reached the range 
of the Ramsauer-Townsend (RT) minimum in the scat- 
tering cross section cr(/id + H2). Then the mean distance 
between collisions increased and the hydrogen layer be- 
came effectively transparent for the d/i (t/i) atoms which 
were then easily emitted from the solid into the adjacent 
vacuum. Such energetic emitted muonic atoms traveled 
through vacuum toward the downstream hydrogen layer 
(DS) and a fraction passed, after possible interactions in 



Au 



Ne 



Au 



18 mm 



initial du 
initial pu 



fast d|r 



slow pu 
\i decay 




FIG. 9: Single target scheme. 



Ne X rays 

FIG. 10: Double target TOF scheme for d/t. Fast and slow d/i 
atoms were emitted forward into the vacuum gap. Although 
not shown in the drawing, muonic atoms were also emitted 
toward the gold foils, which we refer to as backward emission. 
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hydrogen, to the Ne layer and produced x rays as a result 
of the muon transfer to neon and subsequent muonic neon 
deexcitation. Such d/x (tfi) atoms gave a characteristic 
peak in the TOF spectrum at « 1 /is — a time deter- 
mined by the distance between the foils and the position 
in energy of the Ramsauer-Townsend minimum. 

Clearly, not all muons were stopped in the US target 
and thus those which reached the second foil created p/x's 
in the DS protium. The p/x's diffused to neon giving a 
contribution to the time spectrum at early times. Thus 
the resulting Ne time spectrum contained two relatively 
distinct components, one of them connected with the RT 
effect in d/x (t/i) scattering and the other with the dif- 
fusion of p/x atoms in the solid hydrogen. Despite the 
overlap of the two effects in the time spectra one should 
note there is an important advantage of such an experi- 
ment. Since the kinetic energy of the d/x (t/i) is relatively 
high, the delayed RT peak is not sensitive to the state of 
the target material and can be well described using either 
"gas" or "solid" scattering cross sections (see pa H^). 
Due to this, the RT peak can be used as a reference in 
the analysis of the diffusion part where the effects of the 
solid state can be found. 



C. Monte Carlo simulations 

The Monte Carlo code FOW ^3] was used in the plan- 
ning stages of the experiment as well as for the analysis to 
simulate all physical processes occurring after muons pass 
through the entrance window of the apparatus. Muon 
stopping distributions along the beam axis in the dif- 
ferent apparatus components, especially in the hydrogen 
layers of the target, have been taken from a special set 
of measurements |25j and from another Monte Carlo cal- 
culation H3 and used as an input to FOW. The FOW 
code gives the possibility of a full three-dimensional de- 
scription of the target geometry. 

The muonic processes considered are as follows: 

(i) elastic scattering: p/i + p, p/x + d, d/i + p, p/x + t, 
t/i + p, d/i + d, t/i + t. 



(ii) spin-flip: p/i(F) 
t/x(F)^t/x(F'). 



p/i(F'), d/i(F) -> d/i(F'), 



(hi) charge transfer: p/i — > d/i, p/i — > t/i. 

(iv) molecular formation: p/x+p — > pp/i, d/i + d — > dd/i, 
t/i + t — » tt/x, d/i + p — > pd/i, t/i + p — > pt/i. 

Energy-dependent values of the total and differential 
cross sections for the elastic scattering of muonic atoms, 
spin-flip, and charge transfer transitions were used in the 
calculations. For small collision energies (usually below 
0.1 eV) the "solid" double differential cross sections 
were used for the elastic scattering and spin-flip inter- 
actions. At higher energies, where the solid state ef- 
fects become negligible, the total and sin gle differential 
"gas" cross sections from Refs. [9l ITol l4lt l49 f . corrected 



TABLE II: Characteristics of the muonic processes in different 
upstream targets (see column 4 of TableQfor target details) as 
calculated by the Monte Carlo. Muon stops in the targets are 
given in % of muons entering the apparatus (the Monte Carlo 
analog of GMU). The emission yield, molecular formation, 
backward escape, and muon decay are given in % per muon 
stopped. Simulations were performed using the "solid" cross 
sections. 



upstream targets 


DEm 


TEm 


sTEm 


sPP 


PP 


/x-stops 


58.6 


47.4 


32.3 


32.3 


58.6 


pp/x formation 


34.1 


24.3 


4.8 


77.9 


83.8 


pd/x formation 


39.7 










pt/x formation 




48.5 


35.8 






muon decay 


7.1 


6.1 


2.8 


9.5 


10.2 


backward escape 


13.0 


15.2 


31.0 


8.1 


4.1 


forward p/x emission 


1.8 


0.4 


0.4 


4.5 


1.9 


forward d/x emission 


4.3 










forward t/x emission 




5.5 


25.2 







for molecular effects via the Sachs-Teller model, were 
applied. Using "gas" cross sections at higher energies 
saved computer time without incurring any loss in ac- 
curacy. The energy-dependent pd/x, pt/x and dd/x (reso- 
nant and nonresonant) formation rates were taken from 
the Faifman calculations [oOl loll IH^. The pp/x and tt/x 
formation rates were considered as energy independent 
(A PPAI = 3.21 /xs" 1 0, A tt „ = 1.80 /xs" 1 HI). Ther- 
mal motion of the target molecules was also taken into 
account. Tables ITT1 and Hill show the main characteristics 



4000 






b) slowing down 


region 




solid 






gas _r- J V 











p\l energy (eV) 



0.15 0.2 0.25 

p|X energy (eV) 



0.3 0.35 



0.4 
x 10 



FIG. 11: Calculated low energy spectrum (a) of emitted p/x 
atoms for the PP target. The solid line shows the result us- 
ing the "solid" cross sections, whereas the dot-dashed line 
shows the spectrum using the "gas" cross sections. The same 
number of incident muons (2-10 6 ) is used for both cases. Top- 
right picture (b) shows the slowing down energy spectrum in 
a log-binned scale. 
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TABLE III: Characteristics of the muonic processes in down- 
stream protium targets (see column 6 of Table for target 
name references). Muon stops, pfi emission and d/x (t/x) trans- 
mission are given in % of muons passing the entrance window 
of the apparatus. Simulations are performed using the "solid" 
cross sections. 



downstream targets 
(associated US target) 


600 
DEm 


500 
sTEm 


350 
TEm 


300 
DEm 


p-stops 


7.8 


9.3 


2.6 


4.4 


forward p/i emission 


0.6 


0.9 


0.4 


0.7 


dfi transmission 


0.7 






0.9 


t/j transmission 




1.4 


0.9 





calculated for different upstream and downstream targets 
used for emission study. The values are based on at 
least 10 6 simulated muons, so the statistical uncertainty 
is negligible. 

Figure 1111 shows the calculated energy spectrum of 
p/i atoms emitted from the PP target. A strong solid- 
state effect is evident both in yield intensity and spectral 
shapes when comparing the calculations which include 
or disregard the solid hydrogen structure (solid and dot- 
dashed lines, respectively). The details of the spectral 
tail above the Bragg cutoff limit are shown on a log- 
binned scale in the insert. The kinetic energy spectra for 
d/j. and t/i atoms emitted from DEm and TEm targets, 
respectively, are shown in Fig. ^] for the "solid" (solid 
line) and "gas" (dashed line) cross sections. The lack of 
solid state effects is not surprising given the high energy 
of the atoms involved. The muonic atom energy applies 
to atoms which have traveled the 17.9 mm distance be- 
tween the US and DS layers. 

Figure full shows the probability of p/z emission from a 



1 



du 




0.02 



16 20 

energy (eV) 

FIG. 12: Energy spectra of d/j, and tfi atoms emitted from 
DEm and TEm targets, respectively, after flying the US-DS 
distance and before entering the DS target. Solid lines: MC 
with the "solid" cross sections, dashed lines: MC with the 
"gas" cross sections. 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 

(j,-stop position (mm) 

FIG. 13: The probability dependence of p/j, emission from the 
PP target (for a beam momentum p = 26.25 MeV/c) ver- 
sus the initial p/j, position in the target. The simulation was 
performed using both the "solid" cross sections (solid line) 
and "gas" cross sections (dashed line) for the same number of 
incident muons. 



hydrogen target as a function of the initial p^ formation 
position following the muon stop. The example is for 
the PP target and for the muon stopping distribution 
for the beam momentum 26.25 MeV/c. The difference 
between the results of the "solid" cross sections (solid 
line) and "gas" cross sections (dashed line) illustrates 
the strong increase of the mean free path in the final 
stage of the p/i slowing down when the solid state effects 
are considered. One can see that the volume from which 
emitted p/i atoms can originate is much more extended 
in the solid case. 



IV. MEASUREMENTS 

Muonic atom scattering in hydrogen was measured via 
the x-ray time spectra of /«Ne 2p-ls at 207 keV (see 
Fig. I14|l . The time spectra events were selected within 
the energy window 205.6 — 208.3 keV. The background 
time spectra under the 207 keV peak were created us- 
ing time spectra from two neighboring energy windows, 
namely a left background at 203.3 — 205.6 keV and a 
right background at 208.3 - 210.2 keV. Two different 
background evaluation procedures were used. In the first, 
the left and right spectra were added, and then normal- 
ized by the energy window widths and the resulting spec- 
trum subtracted from the time spectrum of the /xNe 2p- 
ls peak. The second relied on a multi-parameter fit of 
the summed left and right backgrounds using two expo- 
nential functions (with lifetimes for muons in gold and 
neon) and the background predicted from the fit func- 
tion was subtracted from the /LtNe 2p-ls peak. Since the 
background accounts for 70 — 85% of the total statistics, 
its removal plays an important role, especially for data 
at early times where the muon prompt capture in the 



10 



10000 




2000 



o 1 — >— ' — 1 

190 200 210 220 230 

X— ray Energy [keV] 

FIG. 14: x-ray energy spectrum for the T3 measurement. 
The fiNe 2p-ls is located at 207 keV, whereas the bigger 
peaks around 215 — 220 keV are the /xAu 6-5 lines. 

neon layer and the gold foils bring a strong contribution. 
Another data cleaning method resulting in better signal 
to background was the requirement that the muon de- 
cay electron be seen after the /iNe 2p-ls x ray, starting 
from a given time delay. Those delayed electrons were 
detected by the scintillators (see Fig. |8J) during a time 
interval 0.2 — 5.2 /is after the /iNe signal. This method, 
called the "del e " criterion, suppressed the background 
by a factor of about 300. However, useful statistics were 
reduced by about a factor ten. 



A. Combined measurements with H/D and H/T 
targets 

For the p/i scattering analysis only the few experiments 
from Table [I] in which the p/i diffusion time spectrum 
have been seen are useful. Similar runs were summed 
whenever possible and Table IIVI give the details. The 
other measurements were nevertheless necessary for the 
determination of detection efficiencies. 

A typical TOF measurement (#1, Table Hvjl . where 
both d/i's from the US layer and p/i's from the DS pro- 
tium layer were detected when they reach the DS neon, is 
presented in Fig. ^| (points with error bars). The a) and 
b) graphs show the time spectra for the full statistics and 
for the events where the del e criterion has been applied 
for background suppression, respectively. That measure- 
ment was performed with a DEm upstream target (with 
0.05% D 2 ) and a 300 Torr (H 2 ) downstream layer. The 
events occurring at early times (t < 600 ns) are due to 
p/x formed directly in the DS hydrogen which then dif- 
fuse to the neon layer. The peak in the TOF spectrum 
corresponds to the delayed d/x atoms which travel the dis- 
tance between the two foils and are not stopped in the 
DS hydrogen due to the RT effect. Also plotted are the 
simulations using the "solid" scattering cross sections as 
well as the result of the calculation when one neglects the 



solid state effects and uses only the "gas" cross sections. 

Another example of a similar TOF measurement (with 
a hydrogen/tritium mixture in the target upstream and 
500 Torr • i protium covering the downstream Ne, #2) is 
shown in Fig. ^] for the full statistics and del e require- 
ment cases. A relatively high muon stopping fraction in 
the downstream target (because the US target was only 
1000 Torr • I and 9.3% of the muons stopped in the DS 
target) gives good statistics for the p/i part in the time 
spectrum. The delayed peak from t/x transfer events lies 
earlier in time than the corresponding deuterium case be- 
cause of the higher t/i energy (see Fig. I12|) . and hence 
the overlap of both diffusion and RT spectra parts is 
fairly strong. The "solid" cross section MC spectrum is 
also presented in the figure (solid line histogram). Dot- 
dashed and dotted lines show the predicted contributions 
from p// and t/i, respectively. 

Figurc lT7l shows the time spectrum of p/x and dfi atoms 
emitted together from a DEm target (#4, Table IIV() . 
This case was analyzed using the different emission time 
dependences since the d/i part of the time spectrum de- 
creases much faster than the pfi emission spectrum; the 
mean diffusion time of d/z in the 2000 Torr ■ I target is 
~ 100 ns, much less than for the pfi emission (~ 300 ns). 



a, 1000 

a 




400 800 1200 1600 




400 800 1200 1600 2000 



time (ns) 

FIG. 15: Experimental time-of-flight spectra (points with 
error bars) for experiment #1 (see Table ITvTl for cases: a) full 
statistics, b) deh criteria. The solid line represents the Monte 
Carlo simulation based on the scattering cross sections when 
solid effects were taken into account, the dotted line is for the 
"gas" cross sections. 
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TABLE IV: Experiments for the p/t + p scattering study. When the hydrogen layer had the same thickness, the different 
measurements were summed to give more statistics. The last column indicates the signals present in the time spectrum. 



No. 


Label 


Hydrogen 


GMU 


Time Spectrum 


# 




[Torr • £] 


xlO 6 




1 


D3+D7 


300 DS 


671 


p/i + delayed d/i from DEm 


2 


T4 


500 DS 


147 


p/i + delayed t/i from sTEm 


3 


T5+T6 


1000 US 


395 


p/i clean spectrum 


4 


D5+D6 


2000 US 


233 


p/i + d/i diffusion from DEm 



The d/i contribution can be removed using the MC sim- 
ulations since d/i emission is independent of solid state 
effects as one can see from Fig.^J The subtraction leaves 
a clean p/i spectrum which can be compared to the MC 
including solid state effects. 

Normalization of the simulated time spectra was done 
using the RT peak which is equally well described by 
both the "solid" and "gas" approaches. The normaliza- 
tion factor n = Y(exp)/Y(MC), where the yields Y{exp) 
and Y(MC) were the total counts in the time interval 
500 — 2500 ns, was calculated from all eight measure- 
ments given in Table [I] which deal with the RT effect. 
The mean weighted value for the germanium detector 
Gl (which was the same in all experimental runs) was 
n = (5.8 ±0.3) x 10~ 4 . This normalization factor, which 
is effectively the detector efficiency, applies to Ne x rays 
detected from the DS neon layer. For experiments #3 



o 

in 



a 
3 
o 
U 



350 - 



300 - 



250 - 



b) del c request 



a) full statistics 




50 - 



-50 



and 4 (see Table IIVII where x rays were detected from 
neon on the US target (shifted by 17.9 mm compared to 
the DS target) the efficiency was « 40% higher (estab- 
lished by comparing the total counts in the prompt Ne 
x-ray peaks in experiments Dl and D5) and the value 
for detector Gl was n = (8.1 ± 1.0) x 10~ 4 . No efficiency 
was determined for the G2 germanium detector, because 
it was changed between the deuterium and tritium mea- 
surements. 



B. Emission of p/i atoms from a layer of solid 
hydrogen 

Although p/i atom emission was a parasitic process in 
the RT experiments, those runs were still useful for the 
analysis of p/i scattering in solid hydrogen. In particular, 
the profound disagreement between experimental data 
and theory using "gas" cross sections and the relative 
correctness of the "solid" cross sections can be seen. 

The observation of enhanced p/i emission from solid 
hydrogen stimulated additional measurements specifi- 
cally intended to study the phenomenon more precisely 
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FIG. 16: Experimental time-of-flight spectra (points with 
error bars) for experiment #2 (see Table HYl for a) full statis- 
tics and b) del e statistics. The solid line is the MC simulation 
based on the scattering cross sections when solid state effects 
were taken into account. Dot-dashed and dotted lines show 
the p/i and t/i contributions, respectively. 



FIG. 17: Experimental time spectrum (points with error bars) 
of p/i and d/t emitted together from DEm upstream target 
(experiment #4, Table ITV1 . Monte Carlo simulation is shown 
as the solid line histogram. Dotted and dashed curves show 
the contributions from p/i and d/t, respectively. Calculations 
based on "solid" cross sections. 
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FIG. 18: Experimental (points with error bars) and MC 
(lines) time spectra of p/i emitted from the 1000 Torr • £ H2 
layer (#3, Table ITVl for the full experimental statistics, a), 
and del e criteria, b). Solid line: calculation with the "solid" 
cross sections; dotted line: calculation for the "gas" cross sec- 
tions. Calculations with "solid" cross sections are normalized 
to the experimental data according to conclusions from the 
TOF measurements. 



with higher statistics. A target similar to the one shown 
in Fig. [§] was made from 1000 Torr ■£ pure protium covered 
with a thin (10 — 20 Torr • £) neon layer. The resulting 
time spectrum of /iNe 2p-ls x rays described the diffu- 
sion of p/i atoms in hydrogen from the moment of the 
muon stop to the moment of emission. 

The measured time spectra (#3, Table HVjl of /iNe 2p- 
ls x rays are shown in Fig. ^| for the full statistics and 
del e criteria. MC spectra also shown on the figures de- 
scribe well the experimental data when the "solid" cross 
sections are used (solid lines). Normalization to the ex- 
perimental data is based on the conclusions from the 

The calculation with 



TOF measurements (see Sec. II V Al 
the "gas" cross sections (dotted lines) gives a suppressed 
yield although the lifetimes representing the diffusion are 
not dramatically different. 



V. ANALYSIS AND DISCUSSION 

As is seen from the comparison between the experi- 
mental results and the simulations performed using both 
the "solid" and "gas" cross sections presented on Fig. ED 



both types of cross section describe the RT part of the 
spectrum equally well. That agreement is due to the fact 
that solid state effects are negligible for the energies of 
the d/i atoms which are responsible for those spectra. A 
similar conclusion can be drawn from the measurements 
with t/i. 

In contrast, there is a big difference in yields as well 
as in time dependence between the calculated p/i diffu- 
sion spectra and measurements where an agreement with 
the experiment is obtained only for the "solid" cross sec- 
tions; the discrepancy seen in Fig. 1151 a) for first points 
from times — 150 ns has an artificial source and will be 
explained later in this section. The "gas" cross sections 
predict a total yield of emitted p/i's two times smaller 
than required for a correct description of the experimen- 
tal results. If one excludes the first channel (Fig. IT5)) . 
which contains events connected with fast p/i atoms from 
the slowing down stage, the "gas" approach gives three 
times less emission. 

The good agreement between experiment and calcu- 
lation based on the "solid" cross sections is also visible 
in Fig. 1161 where the emission spectrum of p/i's from the 
DS target and the transmission spectrum of the delayed 
t/i's from the sTEm US target are presented. The p/i 
emission is enhanced 2.3 times for the total emission and 
and 3.8 times for emission at times t > 30 ns. Accurate 
results, especially from the point of view of the emitted 
p/i diffusion time analysis, are obtained from the exper- 
iment with the sPP target (Fig. I18J) . In that case there 
were no ambiguities between signals from diffusion and 
those from RT events, which was a problem present in the 
measurements with the combined H/D and H/T targets. 
However such combined target measurements were nec- 
essary because they gave the RT peak which was used as 
the reference for the yield normalization of the MC spec- 
tra. The time distribution of the emitted p/i atoms has 
some nontrivial behavior and depends not only on the 
cross sections but also on the target thickness, the muon 
stop distribution in the target, and the initial energy of 
the p/i atoms. In the limit of long times and for a given 
/i-stop distribution, the p/i emission time distribution 
can be well modeled by a one-exponential approxima- 
tion. In that limit, the constant factor in the exponent 
represents the mean time needed for the equilibrated p/i 
atoms to reach the layer boundary (we call that param- 
eter the diffusion time, Td)- Such a Td depends on the 
target thickness but should reach an asymptotic value 
for high thickness simply because the /z-stop distribu- 
tion is effectively exponential decreasing in the thickness 
(c.f. Fig 1 of Ref (43] )• When convoluted with the es- 
cape probability of Fig. and the associated escape 
time from any given depth, this exponential behaviour 
of the stopping yields a constant emission time. In light 
of the remarks presented above, we chose a unique time 
region t > 200 ns for the analysis of the four experimen- 
tal time spectra (see Table Hvjl to determine the diffusion 
time, Td and compare it with the MC simulation. 

Such a choice had the additional advantage that it 
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TABLE V: Diffusion time of p/i atoms, Td in ns, in solid hydro- 
gen layers of different thicknesses (#'s defined in Table ITvt , 



No. Target Experiment MC" 



# 


Torr • i 


full stat. 


del e 


"solid" 


"gas" 


1 


300 


161(32) 


183(35) 


141(4) 


257(26) 


2 


500 


199(40) 


360(165) 


222(7) 


234(14) 


3 


1000 


269(13) 


279(44) 


258(3) 


237(8) 


4 


2000 


279(26) 


253(42) 


254(8) 


230(24) 



"statistical error is used for the fit 



avoided problems with the early parts of the time spec- 
trum. The problem is clearly visible in Fig. 1151 a) for 
times less than 200 ns in a measurement made with a 
thick US layer plus an H2 DS layer. The early time sig- 
nal was only a few percent of the total counts so the 
background subtraction for those time regions was sig- 
nificant, as can be seen from the resulting uncertainties, 
and hence any small irregularities in the thresholds and in 
the detection of secondary muons as well as in the back- 
ground estimation itself can result in an poor estimation 
of the background. It also exists, but is less important, 
in measurements with an H2 US layer, where the per- 
centage of incident muons stopped US was significantly 
higher than in the downstream layer. It is worth noting 
that for the spectra cleaned with the del e condition, and 
hence, where problems with background vanish, the the- 
ory (via Monte Carlo simulations) agrees with the early 
time region, but a more sophisticated analysis was not 
possible because of poor statistics and the complicated 
form of the time spectra. 

The measured and calculated values of Td, fitted using 
a single-exponential distribution, are given in Table 
for both data treatments (i.e., with and without del e ). 
Good agreement between the experimental values of Td 
and the calculations using the "solid" cross sections is ob- 
served for each experiment. The results are also shown 
in Fig. 1191 where the points with error bars represent 
the experimental values of Td from Table [V] The lines 
are the results from the Monte Carlo calculations using 
the "solid" and "gas" cross sections (solid and dotted 
lines, respectively). The MC results from Table IVl were 
obtained strictly for the given experimental conditions 
whereas the calculations represented by the continuous 
lines in Fig. I19l were made assuming pure protium layers 
of increasing thickness and using the same muon stop- 
ping distribution (beam momentum of 26.25 MeV/c) in 
order to show the smooth dependence of Td on the target 
thickness. 

The analysis of the emission yields of the slow p/i atoms 
supports also the use of the "solid" cross sections. This is 
seen by comparing the experimental results with the sim- 
ulations when the "solid" and "gas" cross sections were 
used. The results for the Gl detector are presented in 
Table I VII The comparison was performed for the time 
interval 200 — 600 ns, characteristic for the diffusion pro- 




500 1000 1500 2000 

H 2 target thickness (torr-1) 



FIG. 19: The dependence of the diffusion time Td on H2 
layer thickness. Solid and dashed lines were calculations using 
"solid" and "gas" cross sections, respectively. 



cess, using the efficiencies established from the RT time 
domain (see Section TlV A|) . The agreement between the 
experiments and "solid" cross sections is excellent. 

Despite the general agreement between the experimen- 
tal data and the theoretical description of p/i diffusion in 
fully modeled solid hydrogen, the question of the sensi- 
tivity of the calculated diffusion time and emission yield 
on the cross sections is important. The influence of an 
inaccuracy in the P/1 + D2 (T2) and d/i(t/i) +X (X is any 
hydrogen isotope) cross sections is negligible in the p/i 
diffusion study because of the very small concentration of 
deuterium (tritium) admixtures in the targets. The ques- 
tion of d/x + H 2 , t/i + H 2 cross sections in the Ramsauer- 
Townsend region is discussed elsewhere [2^, |2(J . 

To study the sensitivity of the p/i + H2 scattering sim- 
ulation on the cross sections we performed calculations 
with the "solid" scattering cross sections scaled in the 
low energy region (i.e., collision energy < 0.1 eV). When 
applying a constant scaling factor, between 0.7 and 1.3, 
it only changes the character of the slowing down process 



TABLE VI: Comparison of the calculated and measured (full 
statistics) p/i atom emission yields (in % per GMU) from 
the different solid hydrogen layers for the time interval 200 — 
600 ns. The second column represents the total hydrogen 
thickness in both US and DS. 



No. 


Hydrogen 


Experiment 


MC a 


# 


Thickness 




"solid" 


"gas" 




(Torr • i) 








1 


2000+300 


0.21(2) 


0.225(5) 


0.059(2) 


2 


1000+500 


0.35(6) 


0.345(6) 


0.067(3) 


3 


1000 


0.52(7) 


0.521(7) 


0.313(6) 


4 


2000 


0.38(6) 


0.398(6) 


0.074(3) 



"Only the statistical error is given. 
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by a very small amount at short times. However, it does 
not change the diffusion process in any practical way. 
Variations in the diffusion time Td and the total emission 
yield do not exceed 0.5 — 1% for the mentioned scaling 
range. This is not surprising since the scaling does not 
change the characteristic equilibrium energy of the dif- 
fused p/i which is established by the deceleration and 
acceleration processes (i.e., creation and annihilations of 
phonons, respectively). That energy is close to the Bragg 
cutoff energy Eg, because the coherent phonon creation 
for p/i(F — 0) vanishes below E^ (a weak incoherent scat- 
tering remains) and the inelastic scattering cross sections 
for the acceleration (crf" m ) and deceleration (cr^ ss ) cross 
over near Eb (Fig.EJ). Following these remarks, the value 
of Eg may have an influence on the cold p/i emission (but 
not on the diffusion time). However, the E# value is de- 
fined by the geometrical structure of the hydrogen crystal 
(lattice constant) and thus is known relatively precisely 
(better than 5%). One result of the simulations was that 
such a 5% shift of Es (0.1 meV) gives an increase (or 
decrease) of the p/i emission yield by 1.2%. 

Another important factor which can influence the p/i 
diffusion parameters is the pp/i molecular formation rate, 
which is the most frequent p/i disappearance channel fol- 
lowing muon decay. Both parameters — the diffusion 
time, Td, and the emission yield, Y — are sensitive to 
that process. The value X ppfl = 3.21 ± 0.18 /is -1 
used in our simulations is the more accurate value of the 
two existing measurements performed in solid hydrogen 
to date. Decreasing X pm by one standard deviation re- 
sulted in a 4% increase of the calculated diffusion time 
and a similar increase of the emission yield from a thick 
hydrogen layer. Nevertheless, when one takes into ac- 



count all above systematic errors, the comparison pre- 
sented in Tables and IVII allows us to conclude that 
theory based on the "solid" cross section is still consis- 
tent with the experimental results, contrary to that of 
the "gas" approach. 

In this work we tried to confront the experimental re- 
sults obtained for p/i atom scattering in solid hydrogen 
with the theory of low energy scattering including solid 
state effects. Agreement has been obtained between ex- 
periment and that theory with respect to p/i diffusion. 
In particular, the diffusion time has been found and the 
enhanced yield of p/i emission from the thin solid hydro- 
gen layers has been explained. The results of the study 
proved that the observed solid effects in the scattering 
at low energies (collision energy < 0.1 eV) are correctly 
described by the calculated "solid" scattering cross sec- 
tions. The experimentally observed enhancement in the 
emission of cold p/i atoms creates the possibility of using 
such neutral atomic beams as a tool in further studies of 
muonic processes. 
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